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Unlike typical condensed-matter systems, ultra-cold atoms loaded into optical lattices allow separate control 
of both the particle number and system size. As a consequence, there are two distinct "thermodynamic" limits 
that can be denned for these systems: i) "infinite-volume limit" at constant finite density, and ii) "empty-lattice 
limit" at constant particle number. To probe the difference between these two limits and their crossover, we 
consider a partially occupied lattice and study the transport of non-interacting fermions and fermions interacting 
at the mean-field level into the unoccupied region. In the infinite-volume limit, a finite steady-state current 
emerges. On the other hand, in the empty-lattice limit there is no finite steady-state current. By changing the 
initial filling, we find a smooth crossover between the two limits. Our predictions may be verified using available 
experimental tools and demonstrate a fundamental difference between isolated small systems such as ultra-cold 
atoms and conventional condensed-matter systems. 
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In statistical physics, the thermodynamic limit refers to fix- 
ing the density N p /V - where N p is the particle number and 

V is the volume - and let V — > oo for a continuum system, or 
fixing Np/N - where N is the number of lattice sites - and tak- 
ing N — > oo for a lattice system 1 1 1 . The thermodynamic limit 
is important in the correct description of phase transitions, ex- 
tensiveness of statistical quantities, Bose-Einstein condensa- 
tion, and many other physical phenomena (see, e.g., a recent 
review in Ref. Q). 

In typical condensed-matter systems one cannot arbitrarily 
increase the volume (or number of lattice sites) without simul- 
taneously increasing the particle number. As a consequence, 
one is led to define a unique thermodynamic limit. This is 
not the case for ultra-cold atomic systems, where both N p and 

V (or N) are separately controllable. In this case, we can 
then define two distinct thermodynamic limits: one - that we 
may call "infinite-volume limit" (IVL) - where we fix the den- 
sity while increasing the system size, and the other - we call 
"empty-lattice limit" (ELL) - where we keep the particle num- 
ber constant while increasing the system size. The ELL has 
been less explored in conventional solid-state systems, but it 
is of fundamental importance in statistical physics and may 
have applications in ultra-cold atomic systems [3]|4]. In the 
ELL the Fermi energy is close to the bottom of the energy 
band due to the small number of particles compared to the 
large number of available energy states. Therefore, the physi- 
cal properties in these two limits may be quite different. 

A concrete example of this difference in ultra-cold atoms 
is the strongly-interacting spin-imbalanced (polarized) Fermi 
gas J3] [5] |6j. Its stable structure at low temperatures has 
been proposed, based on many-body theories in the IVL, to 
be either a polarized superfluid or a separated structure of 
a superfluid and a normal phase Q. This agrees with ex- 
perimental results with finite ratios between the population 
of the two hyper- fine states |5J|6). However, when one sin- 



gle fermion is immersed in a cloud of fermions of the oppo- 
site spin, it behaves like a polaron and exhibits very differ- 
ent physics |3|. Due to difficulties in formulating the prob- 
lem with strong interactions, the crossover between the single- 
impurity polaron and the finite-population-ratio superfluid re- 
mains an open question. Finding another case where the tran- 
sition between the two limits can be easily checked, both ex- 
perimentally and theoretically, would be thus desirable. 

In this respect, ultra-cold atomic systems present several 
aspects that differ significantly from those of conventional 
condensed-matter systems, which make them ideal candidates 
to explore such fundamental issues. For one, the interactions 
in ultra-cold atoms can be turned off, thus allowing one to 
experimentally study the role of interactions in physical phe- 
nomena in a controlled way (8). This is of particular impor- 
tance in investigating non-equilibrium properties because, due 
to Coulomb interactions, finding a genuine non-interacting 
system is difficult in conventional electronic systems. In addi- 
tion, it is clearly discussed in Ref. [9] that for a non-interacting 
Bose gas with Bose-Einstein condensation (BEC), the predic- 
tion of the number fluctuations of the condensate as obtained 
from the canonical ensemble does not agree with that from 
the grand-canonical ensemble. This difference is based on the 
fundamental assumption of whether there are exchanges of 
particles with a reservoir, and may be tested in recent exper- 
iments of trapped atomic clouds, which to all practical pur- 
poses can be viewed as isolated systems 1 10 1. 

Here we explore the crossover between the IVL and ELL 
by focusing on one observable: the formation of quasi steady- 
state currents in one-dimensional (ID) atomic gases set out 
of equilibrium in optical lattices. We will focus on non- 
interacting fermions and fermions interacting at the mean- 
field level. It was recently shown that when these ultra-cold 
fermions are loaded into optical lattices and set out of equi- 
librium, a quasi steady-state current (QSSC) - characterized 



by a current plateau as a function of time - may be observed 
ifTTI . In the IVL this QSSC then develops into a true global 
steady state. We can however ask how such a fermionic cur- 
rent evolves from the IVL to the ELL and vice versa. 

For a single fermion in a lattice of size N, we will present 
a general argument which rules out a finite average current 
as N — > oo. For the setup discussed in this paper, and for 
non-interacting fermions, we show explicitly that the current 
decreases as 7V~ 3 / 2 . On the other hand, for N p /N = 0(1) 
on a lattice of reasonable size, a QSSC is readily observable in 
several different setups IfTTI [T21 1. By changing N p and N, the 
QSSC vanishes in a continuous fashion, and no sharp transi- 
tion is found. A finite steady-state current is therefore a many- 
particle phenomenon, not a property of few particles, at least 
for the initial conditions and setup considered in this work. 

We begin with a one-dimensional finite lattice of size N 
and consider single-species fermions first. Due to Pauli ex- 
clusion principle, the s-wave interactions between fermions 
are suppressed so there is virtually no interaction in the sys- 
tem. We model the fermions in the lattice by a tight-binding 
Hamiltonian 

H = -iY J c\c 3 . (1) 

(ij) 

Here (ij) denotes nearest-neighbor pairs, t is the tunneling co- 
efficient, and c\ (Cj) is the creation (annihilation) operator of 
site i. The unit of time is to = h/i. Experimentally, such a ID 
lattice may be realized by inserting an optical barrier |[T3l into 
a ring of optical lattices lfT4ll . The resulting C-shaped lattice 
is geometrically identical to an open ID lattice. Another pos- 
sible technique is the spatial light modulation lfT5l IT6l which 
can produce designed patterns of trap potentials. 

Figure [T] illustrates our proposed experimental setup. Ini- 
tially there is a barrier blocking particles from entering the 
right half of the lattice and the atoms are in the ground state of 
the left half lattice. The system is then driven out of equilib- 
rium by removing the barrier. As the atoms move to the right, 
a current through the middle of the lattice develops. Since 
the system we consider is finite and closed, we employ the 
microcanonical formalism (MCF) ifTTl IT8 1 as implemented in 
Ref. IfTTI 

In the MCF one monitors the evolution of the correlation 
matrix C(t), whose elements are cy(t) = (So\c\(t)cj(t)\So) . 
Here \Sq) denotes the initial quantum state. By using a uni- 
tary transformation Cj — X)fc(^) jkdk> one can rewrite H 
as H = J2 P e pd p d p . Here e p denotes the energy spec- 
trum of H. Explicitly for the Hamiltonian shown in Eq. (TTJ, 

U iP = \f^Ti s ' m (^) and e e p = -2*cos( ^). The index 
p denotes the energy level and should not be confused with 
the momentum. From the equation of motion i(dcj(t)/dt) = 
[cj(t),H] it follows that Cj(t) = T, P ( u )jp d p(°) exp(-ie|t) 
(H = 1 throughout the paper). Then the time evolution is given 
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Figure 1: Schematic plot of the setup. The dashed box denotes a 
barrier at t — so particles can only populate the left half lattice. 
The barrier is then lifted and a current ensues. Here particles may be 
in a superposition of quantum states. 

by 

N 

^(t) = £ (u^Uu^D^ioy^-^. (2) 
p,p'=i 

Here D pp ,(0) = J2i,j( u1 hp( u )p'j c ij( t = °) is the initial 
correlation matrix in the energy basis. 

The current flowing from the left to the right for one species 

is I = -(dN L (t)/dt), where N L (t) = 4 (*)<*(*)• For 

the Hamiltonian considered here, 

/(t) = 2fim{ Cjv/2>JV/2+1 (t)}. (3) 

Experimentally one may prepare several identical setups and 
take density images at different times. The current corre- 
sponds to the rate at which atoms are tranferred to the right. 

We consider the transport of N p single-species fermions 
in a lattice of size N and the filling factor is defined as 
n = N p /N. The initial state has the lowest N p energy states 

of the left half Hamiltonian Hl = — ^Y2i<Uj)<N/2 c l c J oc " 
cupied, while the right half lattice is empty. In the energy ba- 
sis of Hl, the correlation function is D^ p , — 9(N p — p)S pp ', 
where 6(x) = 1 if x > and 0(x) — otherwise. This 
corresponds to a Fermi sea of N p particles on the left half 
lattice with Fermi energy Ep = — 2t cos[ p^fyrr]' ^ ne 
tial correlation function in real space is thus Cy(t = 0) = 
E^=i(Ul) P i(U L ) jp/ D^ if 1 < i,j < (N/2) and zero 

otherwise, where (Ul)]p — J sm ( lv+§) ^ s ^ e un itary 
transformation for Hp. The time-evolved correlation matrix is 
given by Eq. |2| and the corresponding current follows Eq. ((3). 

Figure [2] shows the currents from one particle (N p = 1) 
and 128 particles in a lattice of size N = 1024. Due to fi- 
nite size effects, there is revival behavior when t > (N/2)to 
and we will focus on the physics before the first revival oc- 
curs. In other words, we focus on the finite-time behavior 
(t < (N/2)to) before taking the N — > oo limit. In the oppo- 
site order of limits, t —> oo before N —> oo, the system goes 
through many revivals, and the discussion of the QSSC is not 
meaningful. 
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Figure 2: The currents of non-interacting fermions for (a) N p — 1 
and (b) N p = 128 in a lattice of size N = 1024. 



There is a stark difference in the currents of the two cases. 
A QSSC corresponding to a plateau in the current as a function 
of time can be observed in the N p — 128 case. Similar QSSCs 
have been reported in many different setups ( TO [121 with fill- 
ing of 0(1). In contrast, the current of the single-particle case 
is extremely small and does not exhibit a QSSC. The origin of 
this difference is closely related to how the N —> oo limit is 
approached. As anticipated, when both N p and N are tunable, 
one can i) fix the ratio of n = N p /N as N — > oo, which cor- 
responds to the IVL, or ii) fix N p as N — > oo so that n — » 0, 
which corresponds to the ELL. 

We first prove analytically that for the single-particle case 
the current eventually decays to zero as N —> oo. From 
Eqs Q and {3), for N p = 1 we have 



-2t ^2 sin< i - 11 

p¥=p' 



cos(- 



cos( 



P 7T , 

N + l' 



(4) 



•N + l 

U N/2 , p U N/2+hp ,F(p)F(p'). 

Here F(p) = Y^ll U L , 3 iU jp . In the limit N ->• oo, it be- 
comes F(p) -+ 2\j2 Jq 1 ^ 2 dx[sm(2irx) sm(irxp)] or explic- 
itly, F(p) -> 4y2sin(^)/[(4 - p 2 )n}. Therefore F(p) = 
if p is even, except F(2) = v2/2, and F(p) is finite forp odd. 
In the final expression, all terms with {p = odd,p' = odd} 
cancel each other so only terms with {p = 2,p' = odd} 
and {p 1 — 2,p = odd} contribute and their contributions are 
equal. Thus we consider the contribution from {p = 2,p' = 
odd} and the final result is twice of this contribution. To 
make the analysis explicit, we focus on the long-time limit 
and set t = (N/2)to. One can make other choices as long 
as t/to ~ 0(N/2) and obtain the same conclusion. Due to 
the sin function, the first term oscillates rapidly for large t/to. 
For p = 2, only < p'/(N + 1) < 1/y/N may contribute 
finitely. The current then approaches 
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Therefore as N — > oo (in the ELL), the current I — > 0. 

On the other hand, a finite QSSC in the IVL requires n = 
0(1). This can be seen by observing that for a finite QSSC 
with an averaged value I s to exist during a time interval, say 
from t = (N/A)t to t = (N/2)t , the transmitted particles 



must be AA L - I(t)dt a (A/4)i / s - O(N). 

Since ANl < N p , we conclude that a finite QSSC requires 
n = N p /N = 0(1) as N — > oo, and this is the conventional 
thermodynamic limit in many condensed-matter settings. We 
remark that in this argument one can take any time interval for 
the QSSC as long as it is of order 0(N)to, since one expects 
a QSSC to last for a macroscopic time scale. The linear scal- 
ing of the time duration where a QSSC can be observed with 
respect to the lattice size has been reported in Refs. lfTTlfT2l . 
After presenting the results of the current in the two different 
limits, we now study its behavior in between those two limits. 

From the ELL to the IVL - We first fix the lattice size N and 
study the effect of varying N p . In order not to be distracted by 
small oscillations in the QSSC, we define the averaged current 
as 



(I) = 



1 



(N/l)t 



T + (7V/4)t 



I(t)dt. 



(5) 



To 



In our simulations we chose To = [(A/4) — 3]to to avoid 
the few points around (N/2)to where the current starts to ex- 
hibit boundary effects. When there is a QSSC - like the one 
shown in Fig. |2jb) - (I) gives the magnitude of the QSSC. 
When there is no QSSC, (7) gives an estimate of the aver- 
aged current. Figure |3]shows (I) as a function of N p for fixed 
N = 1024. Due to our selection of the initial ground state, 
1 < N p < 512. 

There is no significant current for N p < 16 so this regime 
is basically in the ELL. A QSSC can be observed for N p > 
64 and we indicate the regime where a QSSC can be found. 
Therefore for N p > 64 the system enters the IVL because n = 
0(1). The vanishing of the QSSC in the transient regime 16 < 
N p < 64 is a slow crossover. Here, the oscillations in the 
current become visible and a plateau of length of 0(N/4)to 
can no longer be identified. We demonstrate this crossover 
by showing the currents (not the averaged ones) for N p — 
16, 32, 64, 128 in Figure[3] The Fermi energy Ep normalized 
to 4t as a function of N p is also shown in the inset of Fig. [3] As 
expected, in the ELL (N p < 16) the Fermi energy is close to 
the bottom of the energy band while in the IVL (N p > 64) the 
Fermi energy is above the bottom of the energy band. Thus 
one can see that a finite QSSC emerges when Ep starts to 
deviate from the bottom of the energy band. 

From the IVL to the ELL - We can, however, follow the 
opposite route, namely fix the particle number N p and vary 
the lattice size N. Figure[4]shows the current at fixed N p = 16 
and different N. Here, we report 32 < N < 8192. We find 
that the current becomes negligible for N > 1024 and the 
system is in the empty-lattice limit. Since the IVL practically 
requires N 3> 1, the system in such a setup cannot reach 
the thermodynamic limit with finite filling. However, a QSSC 
with a plateau of length 0(N/l)to can be observed for 64 < 
N < 128. When N — 32, finite-size effects are strong enough 
that no plateau in the current can be found. Our findings thus 
suggest that although a QSSC should be found for N ^> 1 in 
the most strict sense, it already exists for finite systems with 
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Figure 3: (a) The averaged current (defined in Eq. (|5j) as a function 
of particle number N p in a lattice of size TV = 1024. A QSSC can 
be found when N p > 64. Inset: Normalized Fermi energy Ef/AI as 
a function of N p with TV = 1024. The current as a function of time 

for (b) N p = 16, (c) N p = 32, (d) N p = 64, and (e) N p = 128 is 
shown. 
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Figure 4: (a) The averaged current (defined in Eq. (|5j) as a function 
of N for fixed N p = 16. A finite QSSC can be found when 64 < 
N < 128. Inset: Normalized Fermi energy Ef /At as a function of 
TV with TV P = 16. The current as a function of time for (b) TV = 128, 
(c) TV = 192, (d) TV = 256, and (e) TV = 512 is also shown here. 



reasonable size at filling of order O(l). To clarify how the 
current of the QSSC vanishes as TV increases, the currents for 
N = 128, 192, 256, 512 with TV p = 16 are shown in Figure |4] 
The normalized Fermi energy Ep/(4t) as a function of TV is 
shown in the inset of Fig. |4] Again, one can see that when Ep 
deviates from the bottom of the energy band, a finite QSSC 
starts to form. 

After discussing the IVL and ELL separately, we briefly 
compare the threshold value of the filling factor n for obtain- 
ing a finite current in these two different limits. It is important 
to check that as n increases, a finite QSSC starts to emerge 
when n ~ 1/16 as one can see from Figs. [3] and |4] This is 
based on the premise that the size TV should be large enough 
(TV > 32 in our simulations). The threshold 1/16 corresponds 
to where Ep starts to deviate from the bottom of the energy 
band, and it indicates that the system is large enough to tran- 
sition from the ELL to the IVL. Although the ratio 1/16 may 
be specific to the setup discussed here, a threshold with a fi- 
nite TVp/TV should be a generic feature for a finite QSSC to be 
observable. The finite-size effect introduces another threshold 
for the ELL case when TV becomes too small. Indeed, we ob- 
serve no plateau in the current in the ELL case, Fig. |4] when 
TV < 64 even though n > 1/4. The current shows oscillating 



behavior when TV < 64 due to finite-size effects. To summa- 
rize, there is a common threshold - in correspondence to the 
Fermi energy deviating from the bottom of the energy band - 
at n ~ 1/16 when one tunes either the number of particles or 
the number of lattice sites. However, for the case of decreas- 
ing TV from the ELL limit, there is another threshold when TV 
is too small so there is not enough time for a QSSC to develop 
due to finite-size effects. 

Before closing our discussions of the QSSC in non- 
interacting fermions, we comment on whether it can last long 
enough for experimental observations. The time scale to is 
estimated to be of the order of milliseconds in present exper- 
iments ifTTl . Therefore the plateaus shown in Figs. [3] and [4] 
should be observable since they last for at least tens of to. In 
Ref. IfTTl it has been shown that for non-interacting fermions 
the QSSC is also robust against an additional harmonic poten- 
tial and finite temperature effects, so it may be studied sys- 
tematically in experiments. 

The role of mean-field interactions - Next we address the 
important issue whether the QSSC in the IVL can survive 
when there are interactions among fermions. Here we con- 
sider two-species fermions interacting via weak s-wave scat- 
tering and model the system as the conventional Hubbard 
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model 



N 



Hi, 



= H 



(6) 



Here — c ia ^ la 



Ci a , U is the onsite repulsive coupling con- 
stant, and a is the opposite of a. This model should be appro- 
priate for moderate lattice depth [ 19 1 and the onsite interaction 
may be generated by tuning the system near a Feshbach reso- 
nance [8 1. We assume there are N pa fermions of each species 
on the lattice and define the filling as N pa /N. 

In the presence of interactions, finding the complete set 
of eigenstates for H int for a moderate size of lattices be- 
comes difficult. Therefore, instead of writing down the time- 
evolved correlation matrix in terms of the time-evolved en- 
ergy eigenstates, we evaluate the correlation matrix by solv- 
ing the equations of motion. Explicitly, i(d(d icr Cj a ) /dt) — 
([c] CT , H e ]cj a ) + (cl^Cja, H e \), where [•, •] denotes the com- 



mutator of the corresponding operators. One obtains 



Ma 



dt 



tX a - U {c] s c is c] (J c 1 < J ) + U(c] a c ja c) s c^) 



(7) 



Here X a 

J 



— ( C l+l,<j C j>?) + ( C i-l.cr C J,o-) ( C L C j + l,o-) 

When U/t < 1, we implement the standard Hartree- 

Fock approximation by decomposing (cJg.Cig-cfg.Cjo-) as 
I J ,.. \/,.t 



i c is){c[ a Cj a ). This approximation closes the set of equa- 
tions of motion and we can solve them with the initial con- 

^2p,p' = l(Uint,L)pi(Uint,L)jp' D, 



dition Cij (t 



0) 



int : pp' 



if 



1 < i,j < (N/2) and zero otherwise. Here Df nt ppl is the 
correlation matrix in the energy space of the left-half lattice 
of Hi n t and Ui n t,L is the corresponding unitary transforma- 
tion. For general ratios of N pa /N, it is challenging to find 
(Hj(t — 0). However, in the case where initially there are 
two fermions (of opposite species) per site on the left-half lat- 
tice, Pauli exclusion principle requires D 



(t = 0) = Sn if 1 



L 

int,pp' 



SO 



< hj < (N/2) and zero otherwise. 
We will focus on this case and investigate the effects of in- 
teractions. In our simulations we evolve the two species in a 
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s ) during the evolu- 



symmetric way so that {c' irj Cj 
tion. 

Figure [5] shows the current for U/t — 0, 0.5, 1 with TV = 
128 and N pa = 64. One can see clearly that for each case 
there is a finite QSSC. Therefore the QSSC is robust and 
should be observable in either non-interacting single-species 
fermions or interacting two-species fermions if the interac- 
tion is moderate. The emergence of a QSSC in interacting 
fermions has also been found in different setups using the 
time-depenent density-matrix renormalization group method 
with fixed particle numbers l20l and in interaction-induced 
transport, where time-dependent inhomogeneous interactions 
serve as an internal driving force for the current ETI . 

There are two additional interesting effects shown in Fig. [5] 
The first one is that in the presence of interactions, the plateau 
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Figure 5: (Color online) Currents of one species of two-species 
fermions with interaction strength U/t = 0,0.5,1 (labeled next to 
each curve). Here N = 128 and N pa = 64. 



of the QSSC is smoother and exhibits less noise. The interac- 
tions seem to rectify the QSSC and reduce the modulations on 
top of the plateau. The second is that the height of the plateau 
in the QSSC decreases as the interaction strength increases. 
This implies that even repulsive interactions can slow down 
the transport of fermions propagating from a localized region 
to an empty region. This is consistent with the experimental 
findings of Ref. 12211 . where fermions in a 3D optical lattice 
are initially localized due to a harmonic potential and then re- 
leased from that potential. The experimental data suggest that 
interactions, whether repulsive or attractive, slow down the 
transport of fermions. Here we observe similar phenomena in 
ID optical lattices. 

In summary, we have demonstrated that the current of non- 
interacting ultra-cold fermions on ID optical lattices exhibits 
a smooth crossover as one tunes the ratio of particle number 
and lattice size, N p /N, from the empty -lattice limit to the 
infinite-volume limit and vice versa. The detailed study of 
this crossover is made possible in non-interacting systems us- 
ing the MCF While a finite current is not sustainable in the 
empty-lattice limit, a QSSC showing a plateau in I(t) already 
emerges when the system size is reasonably large with finite 
filling. The threshold value of the filling where the QSSC 
emerges is the same for the two cases, and it corresponds 
to the deviation of Fermi energy from the bottom of the en- 
ergy band. There is another threshold when the lattice size 
is too small and finite-size effects prohibits the formation of 
a QSSC. Moreover, we found that the QSSC also emerges in 
interacting two-component fermions with finite filling. Our 
study advances the understanding of the role of the thermody- 
namic limit in non-equilibrium isolated quantum systems, and 
we hope it will motivate experiments in this direction. 
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